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We trace the origin of the newly determined extragalactic gamma ray background from EGRET 
data to an unresolved population of blazars and neutralino annihilation in cold dark matter halos. 
Using results of high-resolution simulations of cosmic structure formation, we calculate composite 
spectra and compare with the EGRET data. The resulting best-fit value for the neutralino mass is 
m x = 515^75° GeV (systematic errors ~ 30%). 
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The origin of the extragalactic gamma-ray background 
(EGB) has been discussed since the seminal paper on 
gamma ray astrophysics by Morrison in 1958 11. Dif- 
fuse, isotropic gamma-ray background radiation results 
either from the emission of numerous sources too faint to 
be resolved, or from weakly interacting massive particles 
(" WIMPs" ) that have survived as a fossil record of the 
early Universe. The EGRET spark chamber detector on 
board the Compton Gamma Ray Observatory completed 
an all-sky survey above 30 MeV, collecting data from 
1991 until 2000 |2j ■ Subtraction of the foreground plays 
an equally important role in determining the extragalac- 
tic gamma-ray background as in the case of other cosmo- 
logical precision measurements, e.g. the measurements 
of the microwave background and its anisotropics. The 
discovery of a residual galactic gamma ray halo at GeV- 
energies Q prompted improvements of the foreground 
model used in the analysis of the EGRET data. A new 
determination of the intensity of the EGB in the energy 
range of 30 MeV - 50 GeV has been accomplished using 
the numerical code GALPROP for modeling the galac- 
tic gamma-ray foreground, now including an Inverse- 
Compton (IC) component [4|. The EGB spectrum has 
two components: a steep-spectrum power law with in- 
dex a — —2.33 and a strong bump at a few GeV. The 
first analysis of the EGB [5J did not reveal as clearly this 
spectral structure. Guided by the observation that the 
net spectral index of —2.10 ± 0.03 was tantalizingly close 
to the mean spectral index of the resolved extragalac- 
tic EGRET sources (all but Centaurus A and the Large 
Magellanic Cloud are blazars), it was then concluded that 
faint, unresolved blazars were responsible for up to 25% 
or 100% of the background, respectively [a til- Physically 
related sources (such as radio galaxies), large-scale struc- 
tures || or gamma ray bursts |i| could also contribute to 
the EGB. Whatever astrophysical scenario may be con- 
sidered, however, a universal multi-GeV bump resulting 
from the superposition of the spectra of a large, diverse 
population of sources remains suspicious. 
By contrast, the observed energies of the excess bump ap- 
pear naturally in the context of models involving weakly 
interacting, annihilating cold dark matter. MeV-scale 



dark matter particles have recently been discussed as 
a source of the galactic positronium halo [10|. The 
Lee- Weinberg criterion for thermal freeze-out during the 
hot Big Bang, however, renders weakly interacting par- 
ticles with masses much larger than that of the pro- 
ton natural candidates for the cold dark matter [ll|. 
Independently, supersymmetry calls for a new stable 
particle with weak interactions and a mass scale close 
to E F = (1/V2G F ) 1/2 ^ 246 GeV, probably the light- 
est neutralino (x?) 12]. The annihilation of the Majo- 
rana neutralinos in dark matter halos - starting from the 
freeze-out in the hot Big Bang and continuing until the 
present day - produces electromagnetic radiation (along 
with v, p, e) from the decay chains of short-lived heavy 
leptons or quarks. Annihilation lines |13| would only arise 
from the loop-level processes \X 11 an d XX ~~ > Z°7, 
and thus their intensities are generally expected to be 
rather small. The continuum gamma-ray energies are 
kinematically lowered by factors of the order of ten . 
Obviously, the number of dark matter halos must be 
much larger than any possible astrophysical gamma-ray 
source population, and hence the main signature of cos- 
mological neutralino annihilation should actually be a 
rather narrow bump in the EGB at about 10 GeV (lif . 
In this Letter, we show that the observed bump in the 
EGB could well be this signature of dark matter annihila- 
tion, and that there exists an allowed range of neutralino 
candidates in the cosmologically constrained MSSM nat- 
urally explaining this feature when combined with a steep 
astrophysical power-law spectrum component. 



MODELING THE ANNIHILATION 
COMPONENT 

With the differential gamma ray energy distribution 
from jet fragmentation and 7r°-decay df, observed en- 
ergy E and redshift z, the extragalactic gamma ray in- 
tensity due to WIMP annihilation can be written as 
$ 7 (E) = c/47tH x 1/2 (<rv) x ft* M p c 2 rit /m^x 

J^ a lKl±z) 3 /s(E,z)r(z) ttdf E ( 1+ ,)]/£(z) dz , 

[131 Il4 llq where p cr i t is the critical density. Since we 
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consider contributions from annihilations at high red- 
shifts, gamma-ray absorption is included via the attenu- 
ation function k(E, z). For < z < 5 we use the atten- 
uation derived from star formation history [l7|. whereas 
for z > 5 the absorption from interactions with the cos- 
mological relic radiation field [18| is employed. The 
range of integration is limited to < z < 20; gamma 
rays from higher redshifts are negligible. The parameter 
£ (z) is given by £ (z) 2 = M (1 + z) 3 + Ok (1 + z) 2 + Oa- 
In this work, we employ the cosmological "con- 
cordance model" of a flat, dark energy and dark 
matter dominated Universe with the parameters 
(Odm, O m , O k , O a ) = (0.23, 0.27, 0, 0.73). For the dimen- 
sionless Hubble-Parameter h we use the value 0.71 [lflj . 
The annihilation induced intensity scales quadratically 
with the dark matter density and thus strongly de- 
pends on the amount of structure present in the dark 
matter. This dependence is included via the function 
r = l/(/5 2 V) / v p 2 dY ( ~p: mean density over volume V), 
which we use as z-dependently evaluated for cosmologi- 
cal volumes in 20]. Generally speaking, T (z) therefore is 
the "enhancement factor" between a structured universe 
and a completely homogeneous dark matter distribution. 
This enhancement due to structure formation is sensitive 
to the predominant density profile of the dark matter ha- 
los, and therefore subject to some uncertainty. Most high 
resolution N-body simulations yield a universal dark mat- 
ter halo profile p(r) = p s /[(r/r s )^[l + (r/rs) 7 ]^-^], 
where ps and rs denote scale density and radius. Mount- 
ing evidence of the existence of this type of dark matter 
density profile and the validity of the paradigm of hi- 
erarchical structure formation comes from X-ray obser- 
vations of Abell clusters [2lJ and from observations of 
the Lyman-a-forest at high redshifts For our cal- 

culations, we will employ the Navarro, Frenk and White 
(NFW) profile (a = 1, = 3 and 7 = 1) ^ and a 
lower mass cutoff for the halos/subhalos of 10 5 solar 
masses as the baseline case. For the mass-dependent con- 
centration parameter c(M halo , z forma ti on ) = r virlal /rs the 
results presented in [25j are used. This scenario yields 
a present-day enhancement of the flux of 2 x 10 6 , com- 
pared to a completely structureless universe [2(ij . If a 
substantial fraction of the dark matter halos has steeper 
inner slopes, like the Moore et al. profile |26(, the overall 
intensity enhancement might well be a constant factor of 
2-25 larger than assumed here (depending on the inner 
cutoff-radius in case of a singular inner slope oc r or 
steeper) |20). Even steeper inner slopes can arise from 
adiabatic compression by baryons |27|. To account for 
this uncertainty, in this paper we will work with the 
NFW-case of T (0) = 2 x 10 6 , while keeping in mind that 
the intensity could be additionally boosted by a factor 
^ = 0(1. ..10). Substantial clumping of the dark mat- 
ter on mass scales below 10 5 solar masses [28] might 
result in further enhancement of the intensities. We 
compare the EGB intensity due to WIMP annihilations 



TABLE I: Limits of the region of MSSM parameter space that 
have been scanned with DarkSusy for cosmologically interest- 
ing neutralino models not excluded by current accelerator lim- 
its (higgsino mass parameter fi; gaugino mass parameter JX12; 
mass of the cp-odd higgs hia ; ratio of the higgs vacuum expec- 
tation values tan /3; scalar mass parameter ms and trilinear 
soft-breaking parameters for the third generation squarks A t 
and Ab) 
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with EGRET data 0, depending on the neutralino pa- 
rameters (crv) x and m x . The EGRET data points in 
the energy range 50 MeV - 300 MeV are very well de- 
scribed by a power law, presumably due to faint, unre- 
solved active galactic nuclei. The best-fit spectrum is 
7.4 x 10~ 7 x (E/GeV)~ 2 33 phcm" 2 s" 1 sr" 1 GeV^ 1 . A 
steeper spectrum than that of the resolved EGRET 
sources is in fact not unexpected due to the flux- 
spectral-index relation [29j. Adding the annihila- 
tion spectrum to this steep power law, best fit val- 
ues for the cross-section times \t and neutralino 
mass are (<rv) x x = (2.6 ± 0.6) x 10 -24 cm 3 s _1 and 
m x = 515^75° GeV (Fig. la). The inferred neutralino 
mass is independent from the details of cosmic struc- 
ture evolution. To verify that correspondingly high val- 
ues for (cv) x can be obtained within the MSSM frame- 
work while producing cosmologically interesting amounts 
of neutralinos, we use the DarkSusy [3(| numerical rou- 
tines to scan the MSSM parameter space. In Fig. lb, 
we plot valid models that have been found in the re- 
gion of the parameter space described in Table I. In this 
" mA — resonance region" , annihilation resonantly pro- 
ceeds via XX — ► A — »• f f j allowing for a high annihila- 
tion cross-section while still producing the correct relic 
density (^J. There is considerable spread among mod- 
els. In a number of cases, the observed EGB-signature 
can be produced even if ^ is close or equal to unity. 
The models we plot are required to thermally produce 
0.175 > O x h 2 > 0.025. For models producing substan- 
tially less than O x h 2 = 0.1 an additional, non-thermal 
source of neutralinos, e. g. from the decay of heavier relic 
particles, might be considered. For a MSSM-neutralino 
with a mass of 520 GeV, (crv) x = 3.1 x 10 -25 cm 3 s _1 
and a moderate \t of 8 the value of x 2 l v is 0.74, which 
is excellent. The MSSM parameters and resulting EGB 
spectrum for this model are shown in Fig. 2. This neu- 
tralino is gaugino-like (gaugino fraction 0.996) and ther- 
mally produces the correct relic density of O x h 2 0.1. 
In this scenario the mass of the lightest Higgs boson H2 
is 118 GeV. WIMPs with similar mass and cross-section, 
but in other respects different parameters, might, how- 
ever, equally be a possibility. 
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FIG. 1: (a): Results of the x 2 -test of the neutralino annihila- 
tion hypothesis against the measured EGB. The annihilation 
cross-section times ^ is normalized to the NFW-profile case 
(r (0) = 2 x 10 6 ) 

(b): Scatter plot of MSSM neutralinos created by scanning 
the parameter space described in Table I; the rectangle de- 
notes the 520 GeV neutralino further explored in Fig. 2 



COMPARISON WITH ASTROPHYSICAL 
BACKGROUND MODELS 

In order to explain the weak concave behavior of the 
EGB intensity above 1 GeV, as it had emerged from the 
first analysis of the EGRET data [5j, a two component 
nature of the variable blazar spectra was assumed: a 
steep power law as a stationary emission component, and 
a natter power law as a flaring component pj. The pre- 
dicted EGB intensity fit the Sreekumar et al. analysis, 
but does not agree well with the strong bump feature evi- 
dent in the new determination of the EGB (Fig. 2). Con- 
structing a straw persons's model with the same redshift 
evolution, we modify the assumptions of Stecker & Sala- 



mon (SS96) by adopting steeper spectra for the quiescent 
(faint) component, and adding a flatter (flaring) spectral 
component with the hardest spectral index determined 
from EGRET data for a single source, to see how well 
the new result for the EGB can be matched. Evaluating 

$AGN K J n(z)(1 + Z )2[( E(1 + z )/ Eb )-2.33 + 

(E(l + z) /E b ) _1 - 5 ] x k(E, z) with the source density in the 
co-moving frame n(z) cx (1 + z) in the redshift range 
0.03 < z < 1.5, we obtain a coarse, rescaled version of the 
SS96-model, in which the amplitude and break energy Eb 
were chosen to minimize x 2 jv. Details of the luminosity 
evolution are unimportant for this test, and the original 
SS96-curve can be reproduced accurately by choosing ap- 
propriate values for the parameters. The result of fitting 
this straw person's model to the new EGB data is a value 
of 1.05 (Fig. 2). It should be noted that, while the astro- 
physical model can in principle produce an acceptable fit 
to the data, this requires a sharp spectral break at an en- 
ergy Eb of ~ 5 GeV. Blazars, however, have continuously 
varying spectral properties (spectral index, peak ener- 
gies etc.). A sharply bimodal distribution of the gamma- 
ray spectral index - as required here - seems unnatural, 
and the physical origin of the universal crossing energy 
thus mysterious. Moreover, the fraction of sources with 
hard spectra at an energy of 1 GeV at any time would 
have to be about 20% - considerably more than the frac- 
tion of the hard-spectrum sources in the EGRET catalog 
[32J. The blazar model also would imply ~ 1000 sources 
with a > 300 GeV flux of the order of a typical Whipplc- 
source, whereas the steeper power law alone corresponds 
to ~ 40 sources. The first number seems worryingly large 
in view of the ~ 10 confirmed sources, in spite of ex- 
cessive observation campaigns on candidate sources from 
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FIG. 2: Extragalactic gamma-ray background: spectrum as 
determined from EGRET data by Strong et al. (data points) ; 
the upper limit in the (60-100) GeV range is from Sreeku- 
mar et al.; steep power law component (dashed), Stecker & 
Salamon blazar model (dot-dashed), straw person's blazar 
model (dotted line), neutralino annihilation spectrum (orange 
solid line), and combined steep power law plus annihilation 
spectrum (red solid line) 
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radio and X-ray catalogues |33j. Source-intrinsic cutoffs 
well below 100 GeV could, however, remedy this prob- 
lem. The new-generation Cherenkov telescopes H.E.S.S., 
MAGIC, VERITAS, and the GLAST observatory will tell 
the story. 

Discussion: We have arrived at the conclusion that the 
neutralino dark matter scenario is in agreement with the 
observed EGB spectrum. The best-fit value for the neu- 
tralino mass is m x = 515I75 GeV (notable systematic 
errors of ~ 30% can be inferred from the systematic un- 
certainties of the EGRET EGB determination |4J; they 
will be substantially reduced by GLAST). The strength 
of the observed signature can be explained by a combi- 
nation of NFW-type dark matter halo profiles and an 
annihilation cross-section of the order that can be ob- 
tained within the MSSM framework. The rather high 
scale for the mass ladder of the superpartners might ren- 
der direct detection of all but the lightest of these par- 
ticles by the LHC difficult 34]. For the present gener- 
ation of elastic-scatter experiments, the WIMP-nucleon 
cross-sections for the majority of these models are out of 
reach (o" x _ p < 10~ 7 pb), but could be accessible by next- 
generation detectors. For the EGB spectrum to be fitted 
acceptably with blazar models of the Stecker & Salamon 
type, these seem to run into several worrying difficulties: 
A universal spectral crossing energy of a few GeV is re- 
quired by conspiracy, the models predict a higher fraction 
of hard-spectrum sources at GeV energies than observed 
by EGRET, and the models possibly imply a higher num- 
ber of low-redshift blazars above 300 GeV than discov- 
ered with imaging air Cherenkov telescopes. The com- 
bination of astrophysical sources and WIMP annihila- 
tion presented above is an interesting alternative scenario 
that should be probed by the next generation of gamma 
experiments. For the ~ 520 GeV neutralino previously 
discussed, the galactic neutralino population would give 
rise to a faint galactic gamma-ray halo with intensity ~ 
10~ 7 ph cm -2 s _1 sr _1 above 1 GeV. Present data do not 
allow to confirm or rule out such a halo component, but 
for a galactic GeV-component commonly attributed to 
Inverse Compton radiation also neutralino annihilation 
has been proposed as a source |3j, 36] . A robust calcula- 
tion using DarkSusy shows that the corresponding galac- 
tic antiproton flux is not incompatible with the BESS 
measurements 37]. If the astrophysical background can 
be characterised, gamma rays due to neutralino annihi- 
latio n fro m the Galactic Center or external galaxies like 
M87 |38j could possibly be detected with low-threshold 
gamma ray telescopes. For a NFW-profile and the neu- 
tralino presented in Fig. 2, the Galactic Center would 
exhibit a gamma ray luminosity of 5 x 10 35 ergs/s above 
1 GeV (F E >5o GeV = 8 x 10~ n ph cm~ 2 s _1 ) from within 
10 _3 sr. Note that in this scenario the high EGB intensity 
corresponds to a moderate flux from the Galactic Cen- 
ter due to the intricate effects of clumping of the dark 
matter (cf. 39] ) . For the MAGIC telescope at an_£nergy 
threshold of 50 GeV and the halo profile from |4J|, the 
annihilation component in the gamma ray spectrum of 
M87 would be detectable at 5c in 250 hours of observa- 
tion time. - - We thank A. Strong, T. Kneiske, M. Merck 
and R. Rilckl for valuable discussions. Support by BMBF 
(O5CM0MG1) and Helmholtz Gememschaft (VIHKOS) is 
gratefully acknowledged. 



* Electronic address: elsaesser@astro.uni-wuerzburg.de 
[1] P. Morrison, Nuovo Cimento 7, 858 (1958). 
[2] P. L. Nolan et al., IEEE Transactions Nucl. Sci. 39, 993 
(1992). 

[3] D. D. Dixon et al, New Astronomy 3, 539 (1998). 

[4] A. W. Strong, I. V. Moskalenko, and O. Reimer, Astro- 

phys. J. 613, 956 (2004). 
[5] P. Sreekumar et al., Astrophys. J. 494, 523 (1998). 
[6] J. Chiang and R. Mukherjee, Astrophys. J. 496, 752 

(1998). 

[7] F. W. Stecker and M. H. Salamon, Astrophys. J. 464, 
600 (1996). 

[8] S. Gabici and P. Blasi, Astropart. Phys. 19, 679 (2003). 

[9] T. Totani, Astropart. Phys. 11, 451 (1999). 
[10] C. Boehm et al., Phys. Rev. Lett. 92, 101301 (2004). 
[11] H. Pagels and J. R. Primack, Phys. Rev. Lett. 48, 223 
(1982). 

[12] M. K. G. Jungman and K. Griest, Phys. Rept. 267, 195 
(1996). 

[13] L. Bergstrom, J. Edsjo, and P. Ullio, Phys. Rev. Lett. 

87, 251301 (2001). 
[14] P. Ullio, L. Bergstrom, J. Edsjo and C. Lacey, Phys. Rev. 

D 66, 123502 (2002). 
[15] P. Gondolo et al., Proc. of the 4th IDM, York, UK p. 256 

(2003) . 

[16] D. Elsaesser and K. Mannheim, Astropart. Phys. 22, 65 

(2004) . 

[17] T. Kneiske et al., Astron. Astrophys. 413, 807 (2004). 
[18] A. A. Zdziarski and R. Svensson, Astrophys. J. 344, 551 
(1989). 

[19] C. L. Bennett et al., ApJS 148, 1 (2003). 

[20] J. Taylor and J. Silk, MNRAS 339, 505 (2003). 

[21] A. D. Lewis, D. A. Buote, and J. T. Stocke, Astrophys. 
J. 586, 135 (2003). 

[22] R. A. C. Croft et al., Astrophys. J. 520, 1 (1999). 

[23] J. F. Navarro, C. S. Frenk, and S. D. M. White, Astro- 
phys. J. 462, 563 (1996). 

[24] J. F. Navarro, C. S. Frenk, and S. D. M. White, Astro- 
phys. J. 490, 493 (1997). 

[25] V. R. Eke, J. F. Navarro, and M. Steinmetz, Astrophys. 
J. 554, 114 (2001). 

[26] B. Moore et al., Astrophys. J. 499, L5 (1998). 

[27] F. Prada et al., Phys. Rev. Lett. 93, 241301 (2004). 

[28] J. Diemand, B. Moore, and J. Stadel, Nature 433, 389 

(2005) . 

[29] C. von Montigni et al., Astrophys. J. 483, 161 (1997). 
[30] P. Gondolo et al., Proc. of the 4th IDM, York, UK (2003). 
[31] H. Baer et al., JCAP 0408, 005 (2004). 
[32] R. C. Hartman et al., ApJS 123, 79 (1999). 
[33] D. Horan and T. C. Weekes, New Astron. Rev. 48, 527 
(2004). 

[34] M. Battaglia et al., Euro. Phys. J. C 33, 273 (2004). 
[35] R. Aloisio, P. Blasi and A. V. Olinto, JCAP 0405, 007 
(2004). 

[36] W. de Boer et al. (2003), hep-ph/0309029. 

[37] H. Matsunaga et al., Phys. Rev. Lett. 81, 4052 (1998). 

[38] E. A. Baltz et al., Phys. Rev. D 61, 023514 (1999). 

[39] S. Ando (2005), astro-ph/0503006. 

[40] J. C. Tsai, Astrophys. J. 413, L59 (1993). 



